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RELATION OF TURBOJET AND RAMJET COMBUSTION EFFICIENCY TO SECOND-ORDER 
REACTION KINETICS AND FUNDAMENTAL FLAME SPEED 

By J. Howard Childs, Thaine W. Reynolds, and Charles C. Graves 


SUMMARY 

Theoretical studies of the turbojet and ramjet combustion 
process are summarized and the insulting equations are applied 
to experimental data obtained from various combustor tests . The 
theoretical treatment assumes that one step in the over-all chain 
of processes which constitute jet-engine combustion is suffi- 
ciently slow to be the rate-controlling step that determines com- 
bustion efficiency. The parameter P/T/Vy, based on the as- 
sumption that chemical reaction kinetics control combustion 
efficiency , and the parameter P?* 24 Ti/V?* 76 , derived from the 
assumption that the rate of flame spreading controls combustion 
efficiency , are evaluated for various turbojet and ramjet com- 
bustors . In these parameters P t is combustor-inlet static pres- 
sure ) Tj is combustor-inlet static temperature , and V r is com- 
bustor reference velocity based on Pj, T u and maximum com- 
bustor cross-sectional area. 

The parameter P l TJV r has been evaluated for 23 turbojet 
combustors; it provides an adequate correlation of the combustion 
efficiency of most combustors at simulated engine operating 
conditions but proves inadequate in some cases where combustor 
test conditions are varied through wide ranges. The parameter 
po.^Ti/vj?* 76 correlates the data for two small-scale ramjet 
combustors but is generally unsatisfactory for correlating turbojet 
data. Data obtained with one combustor indicate that a shift 
from one rate-controlling step to another occurs as combustor 
pressure is increased; the parameter PjTj/Vj. therefore best 
correlates the data at low pressures , while the parameter 
po. 24 t t/V?* 76 is better at higher pressures. Thus , no single 
correlating parameter can be expected to be adequate for all 
combustors and for the entire range of operating conditions. 

INTRODUCTION 

One of the most serious problems encountered in the 
operation of jet-propelled aircraft is the reduction in com- 
bustion efficiency that occurs at high-altitude flight condi- 
tions. Experimental investigations with both turbojet and 
ramjet combustors have shown that combustion efficiency is 
adversely affected by the high velocities at which these 
combustors are required to operate and by the low pressures 
and low inlet temperatures encountered at high altitudes 
(ref. 1). A theory of the jet-engine combustion process is 
therefore needed in order to explain these effects and to 
indicate the design approaches that are most promising for 
alleviating these adverse effects. This report describes the 
theoretical treatment of the combustion process in turbojet 


and ramjet combustors and the correlation of experimental 
data on combustion efficiency with the theoretical equations. 
The NACA work in this area dates back to 1950 (refs. 2 to 7). 
This work was subject to security classification until 1956, 
at which time much of this material was presented in 
reference 8. 

If the fuel were properly mixed with air and the fuel-air 
mixture were allowed adequate time in the combustor, 
thermodynamic equilibrium would be achieved and the 
combustion efficiency would be 100 percent. The occurrence 
of combustion efficiencies below 100 percent indicates that 
the conversion processes by which the chemical energy of the 
fuel is converted into sensible enthalpy of the exhaust 
products are not rapid enough to proceed to completion dur- 
ing the residence time allowed in high-velocity combustors. 

A theoretical treatment of jet-engine combustion is diffi- 
cult because of the many different conversion processes that 
must occur. The fuel must be vaporized, mixed with air, 
ignited, and oxidized to the final products of combustion. 
These combustion products must then be mixed with dilution 
air to reduce the combustor-exit temperature to the desired 
value. Tile combustion can be visualized as a competition 
between the conversion processes (vaporization, mixing, 
ignition, and oxidation) and the quenching that occurs when 
the reacting mixture is cooled by the dilution air and when 
the reacting mixture contacts the relatively cool walls of 
the combustor. In a ramjet combustor operating at nearly 
stoichiometric fuel-air ratios, dilution-air quenching ceases 
to be important; the efficiency is then determined by the 
extent to which the conversion processes proceed before the 
reacting mixture is swept out of the exhaust nozzle. Because 
of the complicated nature of the over-all process, no exact 
theoretical treatment is possible. 

It is possible, however, to apply theoretical considerations 
in the analysis of the process by use of simplifying assump- 
tions. If the rate of any one of the conversion processes is 
substantially less than the rates of the others, this one process 
will then govern the over-all rate and hence will determine 
the combustion efficiency. The existence of one slow con- 
version step in the over-all chain of processes was assumed. 
Combustion was then treated theoretically, assuming each 
of the various conversion steps to be the over-all rate- 
determining process. The details of the analysis are presented 
herein, and the results are applied to experimental data 
obtained with various turbojet and ramjet combustors. 

1 



2 


REPORT 1334 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


ANALYSIS AND APPLICATION TO COMBUSTOR DATA 

CHEMICAL REACTION KINETICS 

One theoretical analysis (ref. 2) of the turbojet combustion 
process is based on the assumption that the chemical re- 
action (oxidation of the fuel) constitutes the over-all rate- 
determining step in the combustor. This oxidation occurs 
by some chain mechanism. Sometimes, however, the 
kinetics of chain reactions depend entirely on the kinetics of 
a single, slowly occurring reaction within the chain. The 
analysis described herein is therefore based on the kinetics 
of a bimolecular chemical reaction. 

For the bimolecular reaction, A + B-^C + D, in which A 
is the reactant present in smaller quantity, the reaction rate 
is given by the equation 

dt~N A 


(All symbols are defined in appendix A.) Substituting 
appropriate relations from the kinetic theory of gases (ref. 9) 
leads to 

ix i iw 

d t ~p d t~ RHT* [ ) 


Equation (2) applies to any sample of the reacting mixture 
as it passes through the combustor. A detailed derivation 
of equation (2) is given in appendix B and in reference 2. 

At this point in the analysis, two somewhat different 
approaches are taken. These are treated as the nonhomo- 
geneous reactor and the homogeneous reactor, respectively. 

Nonhomogeneous reactor. — Figures 1 and 2 show values 
of fuel-air mixture composition within the burning zone of a 
typical turbojet combustor; these data were obtained by 
a water-cooled sampling probe. The vapor fuel-air ratio 
decreases progressively along the length of the combustor, 
and marked variations in composition occur over the com- 
bustor cross section. The data of figures 1 and 2 indicate 
that homogeneity is not approached. The rate* of mixing 
of fuel and air may nevertheless be quite rapid. Rapid 
fuel-air mixing is assumed in this portion of the paper; 



Figure 1. — Variation of vapor fuel-air ratio with distance from fuel 
nozzle in typical tubular turbojet combustor. Simulated high alti- 
tude; over-all fuel-air ratio, 0.007. 



Figure 2. — Variation of fuel-air ratio with angular position of sampling 
probe 12J4 inches downstream of fuel nozzle in a typical tubular 
turbojet combustor. Fuel, isooctane; over-all fuel-air ratio, 0.014. 


possible effects of slow mixing rates are treated in the section 
on Other Conversion Processes. 

The major assumptions in this analysis may be sum- 
marized as follows: 

(1) Reaction occurs in local zones of stoichiometric fuel- 
air ratio that exist at the interface between the fuel-rich 
and air-rich regions. 

(2) The mixing process whereby fuel and air are introduced 
to the reaction zone is sufficiently rapid that the chemical 
reaction rate controls the over-all process. 

(3) The temperature in the reaction zone is constant along 
the combustor at a value close to the stoichiometric flame 
temperature. 

(4) The mean concentration of each reactant in the 
reaction zone varies along the combustor and is equal to 
the initial concentration minus the fraction X that has 
already reacted at any position in the combustor. 

The foregoing assumptions do not necessarily constitute 
those characterizing the simplest or most plausible physical 
model of the combustion process. Rather, they constitute 
the assumptions necessary to arrive at a final equation that 
fits the experimental data. 

At the outlet of the reaction zone, X is equal to because 
the particular reaction under consideration is the one 
governing the over-all process. Thus, 

X—7} b when l~L and t=L/V av (3) 

Integrating equation (2) with the assumption that the 
temperature in the zone where the reaction A + B— >C+D 
occurs is substantially constant and substituting from 
equation (3) give 


,_JVa 

Ns 71 ” , ^ K :s (<r A +<TB) 2 (N B -N A )ZPe- E/ « T 

1 -% + 2_ VavR^T'l* 


( 4 ) 


The foregoing equation is expressed in terms of the values 
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of P, T, and V av in the reaction zone. In experimental 
investigations, these variables are not measured; the values 
of the combustor-inlet variables P iy T i} and V r are usually 
determined. The values of P, T, and V av must therefore 
be expressed in terms of combustor-inlet conditions. When 
appropriate approximations and substitutions are made 
(appendix C), 

,Na 

N B Vb , „ K 3 (<r A +<r B nN B -N A )e- E ' RT .. y 

m J_ £i/2ys/2 X 


(> 


U 1 


APY 

'Pi) 


^A h G d/A 

P ( ^-T 

t 


nGa\A a ) 


P<T t 


Vr 


( 5 ) 


On the right side of equation (5) the terms are grouped into 
(1) those dependent upon fuel type and fuel-air ratio, (2) 
those dependent upon combustor design, and (3) those de- 
pendent only upon the operating conditions. 

For a given combustor operating with a given fuel at a 
fixed fuel-air ratio and with negligible pressure drop assumed 
across the combustor, equation (5) becomes 


In 


,Na 
1 N B Vb 

1 ~Vb 


-\-K 2 =K^ 


( 6 ) 


The same relation can be expressed as 



( 7 ) 


Figures 3 and 4 show experimental data obtained with an 
annular turbojet combustor (combustor 1) plotted in ac- 
cordance with equations (6) and (7), respectively. A 
straight-line data correlation is obtained in figure 3, as pre- 
dicted by equation (6). Data for a range of fuel-air ratios 
are included in both figures 3 and 4. The theoretical equa- 
tions apply only for a fixed fuel-air ratio; however, for this 
particular combustor the efficiency did not vary appreciably 
with change in fuel-air ratio in the range investigated. 

The parameter PiTJVr has been used extensively at the 
NACA Lewis laboratory to correlate the combustion effi- 



P i r t/Vr, (lb){sec)(°R)/ft 3 

Figure 3. — Experimental data obtained with combustor 1 plotted in 
accordance with equation (6). Fuel, gasoline; variable fuel-air 
ratio. 
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Figure 4. — Experimental data obtained with combustor 1 plotted in 
accordance with equation (7). Fuel, gasoline; variable fuel-air 
ratio. 


ciency of turbojet combustors. Combustion efficiencies 
obtained with 23 experimental and production-model turbo- 
jet combustors have been plotted against PiTilV r . The data 
correlations for 14 turbojet combustors are presented in 
reference 2; the other correlations are unpublished. For 
most combustors the correlation is good; but for others the 
data scatter is great and PiTi/V r cannot be used to predict 
combustion efficiency. For combustors in which combustion 
efficiency varies appreciably with fuel -air ratio, a satis- 
factory correlation can be obtained only by using a different 
correlation curve for each narrow range of fuel-air ratio. 
Illustrative examples of the correlations obtained for two 
other turbojet combustors are shown in figure 5. The data 
for combustor 3 (fig. 5(b)) represent the poorest correlation 
obtained in the many applications of the PiTf/Vr parameter. 

Homogeneous reactor. — When the mixing in the com- 
bustor is assumed to be sufficiently rapid to produce a 
homogeneous mixture throughout the combustor volume, 




(a) Combustor 2. 

(b) Combustor 3. 

Figure 5. — Experimental data for turbojet combustors plotted in 
accordance with equation (7). Fuel, AN-F-28. 
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for given operating conditions the value of X is constant over 
the combustor reaction zone, and it is no longer necessary 
to integrate the equation as in the preceding analytical ap- 
proach. A homogeneous reaction zone was also assumed in 
reference 10 in analyzing the combustion process in a can- 
type combustor having a premixed fuel-air feed and in 
reference 11 in analyzing flames in the wake of bluff bodies. 

For this condition T^Ti+AT s r) b) 


rib 


=(f)G> 

K 3 (ffA + a- B ) 2 NB(l — Vi) (l- 


N A 


~ 1 P 7>-E/R (Ti+AT,r, b ) 

jVb’7 


R^(Ti+AT sVb y' 2 V at 

Rearranging terms yields 


^(r < +Ar 5 ^) 3/2 ^ (T<+A7Wt) 

^(ua+JbJWb 


R l 12 


' ±SL 

S A h C d \A a J 

L t 


PiT \ 


V r 


(8) 


( 9 ) 


Here again the terms on the right side of the equation are 
grouped into (1) those dependent upon fuel type and fuel-air 
ratio, (2) those dependent upon combustor design, and (3) 
those dependent only upon operating conditions. For a 
given combustor, a given fuel, and a single fuel-air ratio, 
equation (9) becomes 


Vb (T i +AZ Vb y' 2 e R(Ti+M _ r ^ P.T, 

~ K * V T 


(1 jVb 5?i> ) 


( 10 ) 


The data for combustor 1 were plotted in accordance 
with equation (10). The fact that the predicted straight- 
line relation was not obtained indicates that the assumption 
of nonhomogeneity (eq. (6)) better fits the experimental 
data for this particular combustor than the assumption of 
homogeneity (eq. (10)). This is as expected, since figures 
1 and 2 indicate that a condition of homogeneity is not ap- 
proached within a typical turbojet combustor. For com- 
bustors having higher rates of mixing in the combustion 
space than those existing in most turbojet and ramjet com- 
bustors, equation (10) may afford a better correlation of the 
data. Equation (10) can also be expressed as 


Vb- 




PiT, 

w : 



(ii) 


and, since T { is small compared with A T s , for a reasonable 
correlation of engineering data for which T t does not vary 
widely, it might be assumed that 




( 7 ) 


This equation is the same as that obtained for the non- 
homogeneous reactor. 

FLAME SPREADING 

Since flame speed is one of the fundamental combustion 
properties that can readily fit into a physical picture of the 
combustion process, a simplified analysis was also made 
assuming combustion efficiency to be controlled by the rate 
of flame spreading. The combustion process was visualized 
as the burning of fuel-air globules of random size and shape, 
which are surrounded by a flame surface and which are con- 
sumed as they pass through the combustor. The flame 
surface advances into the adjacent unburned mixture at a 
rate determined by the physical conditions of the mixture. 
The effect of turbulence on the rate of flame spreading is 
considered solely in terms of its effect on the flame surface 
area. 

The combustion efficiency may be expressed as follows: 


_ p u SU L 

^-pArVr 



The numerator represents the mass of combustible consumed 
by the flame per unit time; the denominator represents the 
mass flow through the combustor per unit time. Hence, the 
ratio of these two quantities must equal the combustion effi- 
ciency. 

Bunsen-burner experiments (refs. 12 to 14) have shown 
that 

^ oc Re°- 2i Pi ' 3 oc S (13) 

u L 


if it is assumed that the effect of turbulence on flame speed 
is due to flame wrinkling (increase in S ). For isooctane- 
air mixtures (ref. 15), 

U l ccT 1a (14) 


The effect of variations in ambient pressure on the laminar 
flame speeds of hydrocarbon fuel-air mixtures has not been 
definitely established. The effect is small and appears to 
vary with fuel type. In general, the relation between 
the flame speed and pressure can be expressed in the form 


U L ozPj 


(15) 


for hydrocarbon fuel-air mixtures. Here q generally has a 
value ranging from —}{ to approximately 0. Reference 16 
concludes that the flame speeds of hydrocarbon fuel-air 
mixtures should be independent of pressure if the tube 
diameter is large enough to prevent surface effects from 
becoming important. In reference 17, the flame speed of 
vaporized isooctane-air mixtures was found to be inversely 
proportional to the cube root of the ambient pressure. 
Reference 18 (p. 177) shows that the exponent q in equation 
(15) depends on the flame speed of the fuel. For fuels having 
flame speeds of the order of that of isooctane, the value of 
q is zero. For propane-air mixtures, reference 19 gives 
q= — 0.3 and reference 20 gives g= — 0.19. A value of 
g— — 0.3 was selected for use in this derivation, since when 
this value is substituted into equation (13) the turbulent 
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flame speed is independent of pressure, a fact which has 
been demonstrated for propane-air mixtures (ref. 13). 
With these substitutions, equation (12) becomes 


r}b=Kt 


P°t'V r 


(16) 


Substituting Re^K^PiVrjTl' 1 into equation (16) gives 



(17) 


The parameter derived from the assumption that flame 
spreading governs the over-all combustion rate did not 
correlate the data for combustor 1, as evidenced by the 
scatter of data points in figure 6. However, most of the data 
points do lie close to a common curve. Thus, the flame- 
speed parameter might possibly correlate data over a part 
of the combustor operating range, and a more critical evalu- 
ation of this parameter appears warranted. 

Data obtained with two 5-inch-diameter ramjet combus- 
tors (ref. 6) are plotted in accordance with equation (17) 
in figure 7. One of these combustors (combustor 4) employed 
a V-gutter flameholder and the other (combustor 5) a con- 
ical flameholder. In both combustors the fuel was injected 
well upstream of the flameholder so that a homogeneous 
fuel-air mixture entered the combustor. The flame-spreading 
parameter of equation (17) gave a reasonably good 



Figure 6. — Experimental data obtained with combustor 1 plotted in 
accordance with equation (17). Fuel, gasoline; variable fuel-air 
ratio. 




(a) Combustor 4. (b) Combustor 5. 

Figure 7. — Experimental data obtained with 5-inch-diameter ramjet 
combustors plotted in accordance with equation (17). Fuels, iso- 
pentane and gasoline. 


correlation of the ramjet combustor data. Data for two 
fuels, gasoline and isopentane, fall on the same curve; this 
is not surprising, since the fundamental flame speeds of these 
fuels do not differ appreciably. 

For variations in fuel type giving sizable variations in U L , 
equation (17) becomes 


Vb : 


U L \ 

~ J \ vr* u L J 


(IS) 


where U L , r is the flame speed of a reference fuel. Gasoline 
was selected as the reference fuel, giving U LfT =lA feet per 
second. Data for combustor 4 operating with 16 different 
fuels (ref. 7) are plotted in accordance with equation (18) in 
figure 8. Only the data for carbon disulfide deviate badly 
from the common curve. 
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OTHER CONVERSION PROCESSES 

In addition to those already discussed, analyses of the jet- 
engine combustion, process were also made with the assump- 
tion that (1) fuel vaporization, (2) fuel-air turbulent mixing, 
and (3) fuel-droplet burning (governed by heat transfer into 
the droplet) were the rate-determining steps in the over-all 
combustion process. 'With each of these three analyses, the 
predicted effects on combustion efficiency of operating vari- 
ables such as combustor-inlet pressure, combustor velocity, 
and combustor-inlet temperature differed markedly from the 
effects observed experimentally; that is, the equations re- 
lating Pi , T u and V r derived in these analyses differed greatly 
from equation (7) and gave no correlation of experimental 
data. The results of one of these analyses have been pub- 
lished in reference 21, which presents the change in burning 
rate of single fuel droplets resulting from a change in oxygen 
concentration of the ambient atmosphere. This change in 
droplet burning rate is much less than that required for this 
phenomenon to account for the observed change in combus- 
tion efficiency with change in inlet oxygen concentration in 
a typical turbojet combustor. 

Since none of these three analyses provided agreement 
between theoretical and experimental data, the tentative 
conclusion may be drawn that neither fuel vaporization, 
fuel-air mixing, nor fuel-droplet burning alone constitutes 
the rate-determining, slowest step in the combustion process. 
The analyses were, however, based on certain simplifying 
assumptions similar to the assumptions included in the 
analyses of the chemical reaction kinetics and flame spread- 
ing, which were previously discussed. Therefore, the pos- 
sibility cannot be overlooked that different assumptions in 
these analyses might result in closer agreement between 
theory and experimental data. 

EXPERIMENTAL VERIFICATION OF THEORY 

VARIATION OF OXYGEN CONCENTRATION 

In order to test the applicability of the equations pre- 
viously derived, a number of special experiments were con- 
ducted. Investigations of mixtures of vapor fuel, oxygen, 
and nitrogen show marked effects of oxygen concentration 
on such fundamental combustion properties as minimum 
spark-ignition energy, quenching distance, and flame speed 
(refs. 15 and 22). Oxygen concentration appreciably affects 
equilibrium flame temperature at stoichiometric or richer 
fuel-air ratios and also affects the concentration of the 
chemical species involved in the combustion; oxygen con- 
centration is therefore a means of varying the kinetics of 
the chemical reactions. 

In combustor tests, variations in oxygen concentrations 
should affect the rate of combustion without appreciably 
changing such factors as inlet velocity, turbulent mixing 
associated with inlet conditions, and fuel -spray character- 
istics. This experimental device should simplify interpreta- 
tion and application of the test data. Accordingly, an 
investigation was conducted (refs. 3 and 4) to determine the 
effect of inlet oxygen concentration on the combustion 
efficiency of a tubular turbojet combustor (combustor 6). 
Combustion-efficiency data were obtained for tbe combustor 
operating both with liquid isooctane and with gaseous 
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Figure 8. — Combustion-efficiency data obtained for V-gutter flame- 
holder in 5-inch-diameter ramjet combustor for 14 pure fuels and 
two blends (ref. 7), plotted in accordance with equation (18). 
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(a) Liquid isooctane fuel. (b) Gaseous propane fuel. 

Figure 9. — Sample of experimental data obtained with combustor 6 
and varying inlet oxygen concentration. Fuel-air ratio, 0.012; 
inlet-air temperature, 40° F. 

propane over a range of combustor -inlet pressures and oxygen 
concentrations. The combustor-inlet temperature and the 
oxygen-nitrogen-mixture mass-flow rate were held constant 
in these tests. 
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A typical set of data obtained for the combustor operating 
with liquid isooctane and gaseous propane is shown in 
figure 9. There are pronounced changes in combustion 
efficiency with variation in both combustor-inlet pressure 
and oxygen concentration. Data of this type were treated 
in terms of both the second-order reaction and the flame- 
speed equations previously derived. 

In the application of equation (5) to the oxygen- 
enrichment data, the burning-zone temperature was arbi- 
trarily taken as the stoichiometric adiabatic equilibrium 
temperature, and the concentration of one of the reactants 
was assumed to be proportional to the oxygen concentration 
a of the combustor-inlet oxygen-nitrogen mixture. Under 
these conditions, the ratio can be assumed constant; 

and, for a given combustor, fuel, and fuel-oxygen-nitrogen 
mixture ratio, equation (5) can be expressed in the form 

Here T eq is the stoichiometric adiabatic equilibrium temper- 
ature, which is a function of oxygen concentration. 

The application of equation (19) to the data of figure 9(a) 
is shown in figure 10. The equilibrium temperatures at the 
various combustor-inlet pressures and oxygen concentrations 
were computed by the methods and charts of reference 23. 
For the combustor data obtained with liquid isooctane, an 
apparent energy of activation E of approximately 37,000 
calories per gram mole satisfactorily correlated the data. 
(In fig. 10 the constant 33,500 is equivalent to E/R; if this 
value is multiplied by i?— 1.987 and divided by 1.8 to 
convert the units, the value of approximately 37,000 is 
obtained for E.) This value of E is in reasonable agreement 
with the apparent energy of activation of 32,000 calories 
per gram mole obtained from adiabatic compression data 
for isooctane-air mixtures (ref. 24). for flame-speed data of 
isooctane-oxygen-nitrogen mixtures, reference 15 found 
that a value of 40,000 gave closer agreement between 
experimental and predicted flame speeds than did 32,000 
calories per gram mole. The effect of fuel-air ratio on the 
form of the correlation curves is shown in figure 11, where 
faired curves are drawn through the correlated data for the 
various fuel-flow rates. 

The data obtained for combustor 6 operating with gaseous 
propane were correlated in a similar fashion. Figure 12 
shows the correlation in terms of equation (19) for the com- 
bustor operating at a fuel-air ratio of 0.012. The value of 
apparent energy of activation E required for best correla- 
tion varied from approximately 27,000 in the low combustion- 
efficiency range to approximately 33,000 calories per gram 
mole in the high combustion-efficiency range. Since the 
scatter of the correlation in the low combustion-efficiency 
range was quite sensitive to the value 7?, a value of E of 
28,000 calories per gram mole was used for figure 12. The 
accuracy of measurements of combustor operating variables 
was not good enough to warrant the determination of a 
more exact value of E between 27,000 and 33,000 calories 
per gram mole. This range of values of E is in approximate 
agreement with that cited for propane in the literature. 
In reference 22, a value of 38,000 calories per gram mole is 



Figure 10. — Experimental data obtained with combustor 6 (fig. 9(a))-, 
plotted in accordance with equation (19). Fuel, isooctane; fuel-air 
ratio, 0.012; inlet-air temperature, 40° F. 



Figure 11. — Effect of fuel-air ratio on correlation of data for combustor 
6. Fuel, isooctane; inlet-air temperature, 40° F. 



Figure 12. — Experimental data obtained with combustor 6 (fig. 9(b)), 
plotted in accordance with equation (19). Fuel, propane; fuel-air 
ratio, 0.012; inlet-air temperature, 40° F. 

given. However, unpublished observations in the applica- 
tion of the Semenov theory to flame-speed data of propane- 
oxygen-nitrogen mixtures indicate that a value of 34,000 
calories per gram mole resulted in better agreement between 
experimental and predicted values. 

For combustor data obtained at constant inlet tempera- 
ture and including variations in the inlet oxygen-nitrogen 
mixture composition, the flame-spreading equation (eq. 
(12)) reduces to 

yo.ir ) ( 2 °) 


439185 — 58 
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Here the small change in Reynolds number resulting from 
the change in viscosity of the inlet mixture with oxygen 
concentration is neglected. 

In reference 15, the laminar flame speeds of isooctane- 
oxygen-nitrogen mixtures at atmospheric pressure and vari- 
ous equivalence ratios were determined over a range of 
initial mixture temperatures and oxygen concentrations. 
The maximum flame speed U L was found to be directly 
proportional to the term a— 12. Similar data were obtained 
for propane-oxygen-nitrogen mixtures (unpublished NACA 
data). For oxygen concentrations of 30 percent by volume 
and below, the flame speed was found to be proportional to 
the term a— 11.5. The constant of 11.5, which represents 
the extrapolated value of a for U L of zero, is in agreement 
with the value of 11.6 cited in reference 25 for the concen- 
tration limit of flammability. Substituting these correlation 
terms for the effect of oxygen concentration on maximum 
flame speed in equation (20) gives the following flame-speed 
equations for constant combustor-inlet temperature: 

For isooctane: 

( 2 .) 

For propane: 

f m- 24 (a- 11.5)1 . 

Vb^J I yoM J 

In these equations the laminar flame speeds of both isooctane 
and propane are assumed to have the same pressure depend- 
ence. The combustion-efficiency data are correlated for 
isooctane and propane in combustor 6 in terms of the pre- 
ceding flame-speed parameters in figures 13 and 14. The 
isooctane data correlate well; however, the propane data 
obtained at the lower combustor pressures fall well below 
the data for higher pressures. 

The reaction-kinetics parameter of equation (19) and the 
flame-spreading parameter of equation (21) both correlate 
the isooctane data obtained with combustor 6. In these 
data only the combustor-inlet pressure and the inlet oxygen 
concentration were varied. Evidently the relative effects of 
P t and a are the same in these equations for the range of 
values investigated in the tests with combustor 6. No clear- 
cut explanation is available for the failure of equation (22) 
to correlate the data obtained with combustor 6 operating 
with propane fuel. However, the combustor was operated 
at lower pressure with propane fuel than with isooctane, and 
the low-pressure data fall farthest from the correlation curve. 

INDEPENDENT VARIATION OF COMBUSTOR-INLET PRESSURE 
AND AIRFLOW RATE OVER WIDE RANGES 

The assumption that chemical reaction kinetics govern 
the over-all combustion rates led to the prediction that 
combustion efficiency could be correlated as a function of 
the parameter PiTJVr , and this parameter is equal to a 
dimensional constant times the parameter Pi/w a . The as- 
sumption that the rate of flame spreading governs the over- 
all combustion rate led to the prediction that combustion 
efficiency should correlate as a function of the parameter 
P°i U Ti/V?' 76 . This parameter is equivalent to a dimensional 
constant times P\'°Ti’ u /w% 76 . The reaction-kinetics analysis 
therefore yields a ratio of exponents of the pressure and 
airflow-rate terms of 2, while the flame-speed analysis yields 



Figure 13. — Experimental data obtained with combustor 6 (fig. 9(a)), 
plotted in accordance with equation (21). Fuel, isooctane; fuel-air 
ratio, 0.012; inlet-air temperature, 40° F. 



Figure 14. — Experimental data obtained with combustor 6 (fig. 9(b) 
plotted in accordance with equation (22). Fuel, propane; fuel-air 
ratio, 0.012; inlet-air temperature, 40° F. 


a ratio of these same exponents of 1.32. Obviously both of 
these parameters cannot adequately correlate the combus- 
tion-efficiency data where the combustor pressure and com- 
bustor airflow rate are varied through wide ranges. How- 
ever, both of these parameters do provide a satisfactory 
correlation of the isooctane data for combustor 6 (figs. 10 
and 13). This is possible only because these data were 
obtained for combustor operation at a single value of airflow 
rate. 

In order to determine conclusively whether the ratio of 
exponents on the pressure and airflow-rate terms should be 
2, 1.32, or perhaps neither of these values, it was necessary 
to conduct a series of combustor tests in which the pressure 



Inlet pressure, in. Hg abs 


Figure 15. — Sample of data obtained with combustor 6 for wide varia- 
tion in inlet pressure and airflow rate. Fuel, propane; fuel-air 
ratio, 0.012; inlet-air temperature, 80° F. 
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(a) Fuel-air ratio, 0.012. (b) Fuel-air ratio, 0.008. (c) Fuel-air ratio, 0.016. 

Figure 16. — Cross plot from faired curves of figure 15 to determine relative effects of inlet pressure and airflow rate on combustion efficiency. 


and airflow rate were independently varied over wide ranges. 
Such an investigation was conducted in turbojet combustor 
6 using gaseous propane fuel and a single-orifice fuel injector 
(ref. 5). A sample of the experimental data is shown in 
figure 15, where combustion efficiency is plotted against 
combustor-inlet pressure with parametric lines of airflow 
rate. The data of figure 15 are for a fuel-air ratio of 0.012; 
similar data were obtained for fuel -air ratios of 0.008 and 
0.016. 

The data of figure 15 are cross-plotted in figure 16(a), 
where combustor airflow rate is plotted against combustor 
pressure with parametric lines of combustion efficiency. 
The slope of the lines in figure 16(a) indicates the ratio of 
exponents on the pressure and airflow rate terms of the 
correlating parameter that best fits the experimental data. 
A transition is evident in the vicinity of a combustor pres- 
sure of 15 inches of mercury. At the very low values of 
combustor pressure and airflow rate, the slope of the curves 
is approximately 2, which* agrees with the slope predicted 
by the parameter PiTJV T . At high values of combustor 
airflow rate and pressure, however, the slope of the curves 



Figure 17. — Experimental data obtained with combustor 6 (fig. 15), 
plotted in accordance with equation (7). Fuel, propane; fuel-air 
ratio, 0.012; inlet-air lemperature, 80° F. 


is of the order of 1.3, which agrees with the value predicted 
by the parameter P°' 24 iy T 7 ?* 76 . 

Similar plots of the data obtained at fuel-air ratios of 
0.008 and 0.016 are shown in figures 16(b) and (c), respec- 
tively. The data of figure 16(b) again indicate a transition 
in the curves, the slope being approximately 2 at low values 
of pressure and airflow rate and approximately 1.3 at the 
higher values. The data for a fuel-air ratio of 0.016 do not 
extend to pressures below about 14 inches of mercury; hence 
the transition in slope of the curves was not obtained in 
figure 16(c). 



Figure 18. — Experimental data obtained with combustor 6 plotted in 
accordance with equation (17). Fuel, propane; fuel-air ratio, 0.012; 
inlet-air temperature, 80° F. 
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The data of figure 16 indicate that, for this particular 
combustor, the parameter PiTJVr will correlate data ob- 
tained at the low values of pressure and airflow rate, while 
the parameter P^^TJVr'™ will correlate the data obtained 
at higher values of pressure and airflow rate. N either param- 
eter, however, adequately correlates data obtained over 
the entire range of pressures and airflow rates. This pre- 
diction is verified by the data of figures 17 and 18, where the 
data of figure 15 are plotted as a function of the parameters 
PiTJVr and P^TJVV\ respectively. 

DISCUSSION 

Although the data of figure 16 are meager and therefore 
insufficient to warrant definite conclusions regarding the 
combustion mechanism, there is indication that a shift from 
one rate-controlling process to some other rate-controlling 
process occurs as combustor operating conditions are varied 
through wide ranges. The data of figure 16, together with 
the analyses presented earlier in this paper, indicate that at 
very low values of pressure in the combustor the chemical 
reaction kinetics are the rate-controlling process. At higher 
pressures, however, the rate of flame spreading appears to be 
the rate-determining process. This trend of events is at 
least in qualitative accord with expectations, for the chemical 
reaction kinetics are a function of the square of combustor 
pressure and would therefore be expected to be very slow at 
low-pressure conditions and increase rapidly with increase in 
pressure until some pressure is reached above which kinetics 
no longer constitute the slowest step in the over- all chain of 
processes. More extensive investigations with other com- 
bustors would be required, however, before any definite 
conclusions could be made in this regard. 

Other investigators (e. g., refs. 26 and 27) have obtained 
correlating parameters different from those presented herein. 
In their correlating parameters the ratio of exponents on the 
pressure and airflow-rate terms lies between the values of 1.3 
and 2 obtained for the two parameters in this paper. This 
can be explained if it is assumed that reaction kinetics are 



Figure 19. — Performance curves for three hypothetical combustors. 
For combustor 7, P\jw a should correlate v data. For combustor 9, 
Pf- 24 /w®* 76 should correlate v data. For combustor 8, Pi/w af where 
1.3 < Cw < C2, should provide best correlation of V data obtainable with 
a single parameter; however, no single parameter is really adequate 
for this combustor. 


rate-controlling at low pressures and flame spreading is 
controlling at higher pressures, as follows: Consider figure 19, 
which shows operating curves for three hypothetical com- 
bustors 7, 8, and 9. Combustor 7 gives substantially 100- 
percent combustion efficiency throughout the high-pressure 
range where flame spreading controls the combustion rate; 
its efficiency drops below 100 percent only at the low pres- 
sures where reaction kinetics control the over-all rate of 
combustion. Hence, the data for this combustor would be 
expected to correlate with the reaction-kinetics parameter 
PiTJVr or pyw a . Combustor 9, on the other hand, will not 
operate at low pressures, and its performance should there- 
fore correlate with the flame-spreading parameter P\ U TJ 
Vr' 7 \ which is equivalent to Pi°n u /w° a 76 , or P\' Z2 /w a for 
tests where T t is constant. Combustor 8 gives low efficiencies 
at pressure levels where both flame spreading and reaction 
kinetics are controlling. Its performance data cannot be 
correlated with either of the preceding parameters by itself. 
For a perfect correlation the kinetics parameter should be 
used for the lower operating pressures and the flame-spreading 
parameter for the higher pressures. If all the data are 
plotted against a single parameter, however, then the best 
correlation should be obtained with some parameter Pifw a , 
where 1.32 <Cn<j2. Thus, this argument can account for 
the fact that each investigator (e. g., refs. 26 and 27) seems 
to obtain a different correlating parameter for the particular 
combustor investigated. 

The hypothesis that reaction kinetics are rate-controlling 
at low pressures and flame spreading is controlling at higher 
pressures also gives some insight into the reasons why the 
reaction-kinetics parameter better correlates most turbojet 
combustor data while the flame-spreading parameter is better 
for the small-scale ramjets of figure 7. These ramjets could 
not operate at combustor pressures below 15 inches of 
mercury, while many of the turbojet combustors have been 
operated to much lower pressures. The story is not quite 
this straightforward, however, for in some turbojet com- 
bustors which gave low efficiencies at pressures above 15 
inches of mercury the reaction-kinetics parameter has given 
better data correlations than the flame-spreading parameter. 
Perhaps the pressure level at which the transition occurs from 
one rate-controlling process to the other is appreciably dif- 
ferent in different combustors. 

CONCLUDING REMARKS 

The parameter PiTi/V r , which was derived from the as- 
sumption that chemical reaction kinetics control combustor 
performance, provides an adequate correlation of the 
combustion-efficiency data of most turbojet combustors at 
simulated engine operating conditions. Many of these same 
data are not correlated by the parameter P?‘ 24 7yF? >76 , de- 
rived from the assumption that the rate of flame spreading 
controls combustor performance. The flame-spreading pa- 
rameter does, however, correlate data obtained with two small 
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ramjet combustors. Data obtained with one turbojet com- 
bustor indicated that a shift from one rate-controlling step 
to another occurs as combustor pressure is increased, and the 
parameter PiTJVr therefore best correlates the data at low 
pressures while the parameter P?* 24 T t /Vr' 76 is better at higher 
pressures. No single correlating parameter can therefore be 
expected to be adequate for all combustors and for the entire 
range of operating conditions. 

Combustor performance data are not available for wide 
ranges of inlet-air temperature. The maximum variation in 
Tt was 50 percent for the test data used to evaluate the 
correlating parameter. Consequently, some doubt exists 
whether the parameters derived herein would be adequate if 
larger T t variations were encountered. 

Because of the complexity of the over-all process of turbine- 
engine combustion, it appears unlikely that this simplified 
theoretical treatment (which assumes one step in the conver- 


sion processes to be rate-controlling) would be adequate for 
all combustors. Previous investigations have shown that 
the parameter PiTJVr is not adequate for correlating the 
data for certain turbojet combustors, even for the limited 
range of conditions required to simulate flight operation. 
In addition, experimental data are available that indicate 
fuel vaporization and fuel-air mixing play an important role 
in some combustors (unpublished NACA data). Conse- 
quently, the equations derived in this paper must be applied 
judiciously and with due regard for their limitation. 

The practical use of the parameter PiTi/V r for predicting 
the performance of new combustor designs is discussed in 
reference 28. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 5, 1957 
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APPENDIX A 


SYMBOLS 


average cross-sectional area of combustor reac- 
tion zone, sq ft 
liner air-entry hole area, sq ft 
hole area in primary zone, sq ft 

total hole area in liner, sq ft 

maximum cross-sectional area of combustor flow 
passage, sq ft 

discharge coefficient of liner holes 
energy of activation for reaction, ft-lb/lb 

constants 

length of reaction zone, ft 

distance from upstream end of reaction zone to 
plane of any point within reaction zone, ft 
number of molecules of reactants A and B per 
pound of original mixture, lb” 1 
static pressure in combustor reaction zone, lb/sq 
ft abs 

combustor-inlet static pressure, lb/sq ft abs . 
static-pressure drop across combustor, lb/sq ft 
abs 

gas constant, ft-lb/(lb)(°R) 

Reynolds number 

total flame surface area in combustor, sq ft 
absolute static temperature in combustor reac- 
tion zone, °R 

equilibrium flame temperature, °R 
combustor-inlet absolute static temperature, °R 


A T s 
t 

U L 

U T 

V av 

V r 

w 

dW/dt 

W a 

X 

z 


Vb 

p 

Pu 


<P 


temperature rise for complete combustion of 
stoichiometric mixture, °F 
time, sec 

laminar flame speed, ft/sec 
turbulent flame speed, ft/sec 
velocity based on density, mass flow, and average 
cross-sectional area of reaction zone, ft/sec 
reference velocity based on combustor-inlet den- 
sity, total airflow, and maximum cross- 
sectional area of combustor flow passage, ft/sec 
weight of mixture reacted at any time t , lb/cu ft 
reaction rate, lb mixture reacting/(sec)(cu ft) 
combustor airflow rate, lb/sec 
fraction of N reacted at any time l or at any loca- 
tion within reaction zone 
number of molecular collisions of type involved 
in reaction, per unit time, ft” 3 sec _1 
volume-percent oxygen concentration at com- 
bustor inlet 

probability factor, fraction of collisions involving 
sufficient energy for reaction to occur which 
actually result in reaction 
combustion efficiency 
density, lb/cu ft 

density of unburned gas, lb/cu ft 
effective molecular diameter of reactants A and 
B, ft 

fraction of collisions involving sufficient energy 
for reaction to occur 
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APPENDIX B 


SIMPLIFIED DERIVATION OF KINETICS EQUATIONS 


For the bimolecular reaction, A + B— >C + D, in which A 
is the reactant present in smaller quantity, the reaction rate 
is given by the equation 

d W_Zpp n 

dt~N A 1 ’ 

The probability factor (3 will be assumed to be constant, 
which is conventional. 

A molecular collision of the required type (A + B) occurs 
when the distance between the centers of the molecules is 
(c- a + 0 - b )/ 2 . The average A molecule will therefore sweep 
out in unit time a cylindrical space of length u having 
an effective cross section of (Tr/4)(a A -\-a B ) 2 . The number of 
collisions of the average A molecule with B molecules per 
unit time will be (w/4)(o- A -[-<T B ) 2 un B . Each of the other A 
molecules is also colliding with B molecules; the total num- 
ber of A + B collisions per unit time per unit volume is 
therefore 

Z=“ (a A -\- a B) 2 un A riB 


Substituting V K 7 RT for u, pN A ( 1— X)for n A , and pNs 
^1— for ns gives 

Z=K 3 (c A +<r B y(RT)X(N A NB)(l-X) (l-^ X) p 2 
When PJRT is substituted for p, 

s:,(t,+« b )w»n b (i-X) ( i x ) /" 

z- w* v - 1 (B2) 

From the Maxwell-Boltzmann distribution law (ref. 9), 


p—e~ E/RT 


(B3) 


Equations (B3) and (B4) are substituted in (1) to obtain 


- E/RT 


d t 




CB4) 


Certain approximations are involved in this derivation; 
and, when corrections are made for the distribution of 
molecular velocities and other factors, the equation becomes 

Z=K X (<r A + <T b) 2 un A n B (B 1 ) 


Now 

&X 1 dW *.<«+«>w.(i-x)(i 
d t P d t R*T* 


APPENDIX C 

KINETICS EQUATION EXPRESSED IN TERMS OF MEASURED QUANTITIES 


Equations (2) and (4) are expressed in terms of quantities 
P, T, and V av in the combustion zone. In most experi- 
mental investigations these quantities are not measured. It 
is therefore desirable to express these equations in terms of 
the measured quantities P if T t) and V r . This can be done 
as follows: 

p = p < O-f) < oi) 

where AP/P*«0.05 for most turbojet combustors; 



where A r /A a & 1.7 for most turbojet combustors. Substi- 
tuting these expressions into equation (4) gives 


In 


1 


__N A 
Nb 176 


1 — Vt 
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